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The Crab pulsar is the only astronomical pulsed source detected above 100 GeV. The emission mechanism of 
very high energy gamma-ray pulsation is not yet fully understood, although several theoretical models have 
been proposed. In order to test the new models, we measured the light curve and the spectra of the Crab pulsar 
with high precision by means of deep observations. We analyzed 135 hours of selected MAGIC data taken 
between 2009 and 2013 in stereoscopic mode. In order to discuss the spectral shape in connection with lower 
energies, 4.6 years of Fermi-LAT data were also analyzed. The known two pulses per period were detected with 
a significance of 8.0 cr and 12.6 a. In addition, significant bridge emission was found between the two pulses 
with 6.2 (T. This emission can not be explained with the existing theories. These data can be used for testing 
new theoretical models. 


1. Introduction 


The Crab pulsar and the surrounding Crab nebula 
are the remnant of the supernova of AD 1054. Both 
the pulsar and the nebula are well studied in a very 
wide energy range starting from radio (10“^ eV) to 
VHE energies ( up to tens of TeV). It is one of the 
youngest pulsars known and its spin down luminos¬ 
ity (4.6 X 10^® erg/s) is the highest among Galactic 
neutron stars. To date, this pulsar is the only one for 
which pulsed emission has been detected above 100 
GeV. 


Gamma-ray pulsation from the Crab pulsar up 
to ^ 10 GeV h ad been known since the 1990s 
Nolan et al.l Il993 |. In 2008, pulsations were de¬ 


tected by the MAGIC te lescope at energies above 
25 GeV [Aliu et abl 12008 1. This result suggested 
that the emission originates in the outer magneto¬ 
sphere. The simplest curvature radiation scenario 
in the outer magnetosphere predicts an exponen¬ 
tial cutoff in the energy spectrum at GeV energies 


|e.g., Mu slimov and Hardingjl2004 iTakata et al.ll2005 
iTang et 'al.fl200^ . Fermi-LAT observations from 100 
MeV to a few tens of GeV, which started in Au¬ 
gust 2008 . showed a clear break in the spectrum at 


6 GeV 


Abdo et ^l2010j| supporting this scenario. 


A few years later, h owever, MAGIC and VERITAS 
Aleksib et al.l l201lL l2012aL lAliu et aT] 12011 found 


that the energy spectrum of the Crab pulsar extends 
up to 400 GeV following a power law. The emission 
above 100 GeV is difficult to explain only with the cur¬ 
vature radiation, and additional or different emission 
mechanisms are required. Several new models were re¬ 
cently proposed t o explain the energy spectrum of the 
Grab pulsar [e.g.. lAleksic et al.ll201lLlAharonian et al.l 


l2012| . 

Here we present new results from the continuing 
monitoring of the Crab pulsar with the MAGIC tele¬ 
scopes that will help to constrain any model for the 
emission. In order to discuss the Crab pulsar spec¬ 
tra at energies lower than those accessible to MAGIC, 
Fermi-hAT data were also analyzed. 


2. Instruments, data sets, and analysis 
methods 

2.1. The MAGIC Telescopes 


The MAGIC telescopes are two Imaging Atmo¬ 
spheric Cherenkov Telescopes located on the island 
of La Palma (Spain) at 2200 m above sea level. Both 
telescopes consist of a 17 m diameter reflector and a 
fast imaging camera with a field of view of 3.5°. The 
trigger threshold for regular observations at zenith an¬ 
gles below 35° is around 50 GeV and the sensitivity 
above 290 GeV (in 50 h) is 0.8% of the Grab neb- 
ula flux with an ang ular resolution better than 0.07° 
[Aleksic et al.ll^I2i| . 

For this study we used 135 hours of data taken at 
zenith angles below 35° during optimal technical and 
weather conditions between September 2009 and April 
2013. Standa. r d MA GIC analysis, as describ ed in 
iMoraleio et al.l |2009l| and lAleksic et al.l |20I2bl| , was 
applied to the data. The conversion from event arrival 
times to pulsar rotation al phases used Tempo2 soft¬ 
ware iHobbs et al. l2006| and a dedicated package in¬ 
side MARS LopezIbO^ . The spin parameters of the 
Crab pulsar were taken from the monthly reports of 
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the Jodrell Bank Radio telescope^ [Lvne et ^Il993 |. 3. ReSUltS 


2.2. Fermi-LAJ 


The Large Area Telescope (LAT) is a pair conver¬ 
sion gamma-ray detecto r on board the Fermi satel¬ 
lite [Atwood et al.l 1200^ . It can detect high-energy 
gamma rays from 20 MeV to more than 300 GeV. It 
has been operational since August 2008 and all the 
collected data are publicly available. In this work, we 
have used 5.5 years of Pass 7 reprocessed data^ from 
2008 August 4 to 2014 January 31. The region of in¬ 
terest was chosen to be 30° around the Crab pulsar. 

Along with the public data, the LAT team pro¬ 
vides the corresponding analysis software and instru¬ 
ment response functions (IRF) designed for the anal¬ 
ysis of that particular dataset. We have used the ver¬ 
sion v9r32p5 of the Fermz-LAT ScienceTools^ and the 
P7REP_SOURCE_V15 IRF. From the downloaded 
data we have discarded events taken at zenith an¬ 
gles above 100° to reduce the contamination of albedo 
gamma rays coming from the Earth’s limb. To com¬ 
pute the pulse phase, we used the same spin parame¬ 
ters as for the MAGIC analysis. The obtained fluxes 
were computed by maximizing the likelihood of a 
given source model using the gtlike tools. The binned 
likelihood method was adopted and a 40° square area 
with 0.2° bin width was used for the likelihood max¬ 
imization. Apart from the Galactic (galJem_v05.fits) 
and extragalactic (iso_source_v05.txt) diffuse emis¬ 
sion, we considered as background sources for the like¬ 
lihood fits all sources listed i n the second LAT source 
catalogue Nolan et'^l2012 |. The data taken during 
the periods when the Crab nebula was flaring were 
not excluded from the analysis. These flares should 
not have any impact on the pulsed emission results 
because it is known that the pulsation componen t did 
not change during the flares Buehler et al.ll201^ , and 
the average nebula flux including flare periods was 
subtracted when the pulsar signal was determined. 
Regarding the reported Fermi-LAT spectrum from the 
Crab nebula, the six Crab flares that lasted a few days 
might be responsible for a few percent of the photons 
below 1 GeV in the overall 5.5 year dataset. Given 
that the effect is expected to be small, and that this 
paper focusses on the emission from the pulsar, we did 
not correct for this effect. 


3.1. Light curve above 50 GeV 

Figure [T] shows the light curves of the Crab pulsar 
measured by MAGIC. Two peaks are clearly vis ible. 
Following our previous study [Aleksic et al.ll201^ , we 
define phase ranges for the two peaks as PIm (phase 
—0.017 to 0.026) and P2 m (0.377 to 0.422). The back¬ 
ground level (hadrons and continuum gamma rays) is 
estimated using the phase range between 0.52 and 0.87 
and it is then subtracted from the histograms"^. The 
number of excess events in PIm between 50 GeV and 
400 GeV is 930 ± 120 (8.0 a) and in P2 m is 1510 ± 120 
(12.6 cr). 

In addition to the two main peaks, significant emis¬ 
sion between them is also visible. The region between 
the peaks is generally called the Bridge. Defining 
the Bridge region as the gap between PIm and P2m, 
namely, between 0.026 and 0.377 (hereafter BridgeM), 
we obtain an excess of 2720 ± 440 (6.2 a) events in 
this region. Adopting the definition used at lower en- 
ergies for the Brid ge as the region 0.14 — 0.25 from 
iFierro et al.l |1998 | (hereafter Bridges), then the num¬ 
ber of excess events is 880 ± 200 (4.4 a). This excess 
increases to 1940 ± 370 (5.2 cr) if we extend Bridges 
with the so-called trailing wing of PI and the lead¬ 
ing wing of P2, namel y to the interval of 0.04 — 0.32 
[see lFierro et al.iri998|| . It should be noted that this 
detection co nfirms the hint of br idge emission already 
reported in [Aleksic et al.l[20T^ |. 


3.2. Comparison with iower energies 

Figure [5] shows the light curves at optical. X-ray, 
and gamma-ray energies obtained with various instru¬ 
ments, together with the 50—400 GeV light curve from 
the bottom panel of Fig. [T] The background was sub¬ 
tracted in the same way as the MAGIC light curves 
(see Sect. 13.11) . The intensity and morphology of the 
bridge emission varies considerably with energy. It is 
very weak at optical wavelengths and in the 100 — 300 
MeV range, while there is an appreciable difference at 
X-rays and soft gamma rays. At the energies covered 
by MAGIC, the peaks become much sharper and a 
prominent bridge emission appears. 

It is known that the flux ratio between the two 
peaks strongly depends on energy, as does the ra- 
tio between the fir st peak and the bridge [see, e.g., 
iKuiner et al.ll200l |. Fig. [3] shows the flux ratio be¬ 
tween P2m and PIm and that between BridgeE and 
PIm as a function of energy from optical 2 eV) 


^http: //www .jb.man.ac.uk/~pulsar/crab.html 
^http: //fermi.gsfc.nasa.gov/ssc/data/analysis/ 
documentation/Pass7REP_usage.html 

^ http: / / fermi.gsfc.nasa.gov/ssc / data/analysis / scitools / 
overview.html 


^ An estimation of the background using the off-peak interval 
from the LAT Second Pulsar Catalog, namely the phase range 
between 0.61 and 0.89, lead to very similar results. 
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Figure 1: Light curves of the Crab pulsar obtained by MAGIC from 50 GeV to 100 GeV (top), from 100 GeV to 
400 GeV (middle), and for the full analyzed energy range (bottom). The bin widths around the peaks are 4 times 
smaller (0.005) than the rest (0.02) in order to highlight the sharpness of the peaks. 


to 400 GeV. Steady emission was subtracted before 
the ratios were computed. The ratios P2m/P1m and 
BridgeE/PlM behave similarly. These ratios increase 
with energy up to 1 MeV, decrease up to 100 MeV, 
and increase again from that energy on. At 50 — 400 
GeV, the ratios basically follow the trend seen at lower 
energies. 


whereas the systematic error of the flux normalization 
is estimated to be 18%. We estimate the overall sys¬ 
tematic uncertainty uncertainty on the spectral slope 
to be 0.3. 


4. Discussion 


3.3. Spectrai energy distribution 

The spectral energy distributions (SEDs) of the 
PIm, P2mj BridgeM, and Bridges between 100 MeV 
and 400 GeV are shown in Fig. 01 together with the 
Grab nebula SED obtained with a subset of the data 
used for the pulsar analysis. The SEDs were calcu¬ 
lated using Fermi-LAT data below 50 GeV (below 
200 GeV for the nebula), and MAGIC data above 50 
GeV. The nebula SED is connected smoothly between 
the two instruments. The Fermf-LAT data were fit 
with a power law with an exponential cutoff, while the 
MAGIC data were fit with a simple power-law func¬ 
tion. The obtained fit parameters are summarized in 
Table H The power-law indices between 50 GeV and 
400 GeV are about 3 and no significant difference is 
seen between different pulse phases. The uncertainty 
in the absolute energy scale is estimated as 17%, 


In summary, the Crab pulsar above 50 GeV ex¬ 
hibits a light curve with a significant bridge emission 
between two sharp peaks (Fig. 1). The flux ratios 
P2m/P1m and BridgeE/PlM increase with increasing 
photon energy between 100 MeV and 400 GeV (Figs. 2 
and 3). Between 30 GeV and 400 GeV, the fluence in 
the bridge phase is comparable to that in the PI phase 
(Fig. 4). The SEDs in the 50 — 400 GeV range could 
be fit with power-law functions for the three phases. 

There are several models which can ex- 


such 


as Aleksic et al.l 

20111. Aharonian et al. 

BednarekI l2012l. 

Arka and DubmJ (2013l|. 


can explain the VHE pulse profile and the spetrum 
consistently. Further theoreticaly studies and deeper 
observations of the Crab and other gamma-ray 
pulsars are needed to understand the VHE emission 
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Table I Spectral Parameters 


phase 


■n a 

i 1 

[GeV] 

-^100^ [10“^^TeV“^cm“2s“l] 

F f) 

i 2 

PIm 

8.87 ±0.14 

1.88 ±0.01 

3.74 ±0.15 

4.18 ±0.59 

3.25 ±0.39 

P2m 

3.14 ±0.07 

1.97 ±0.01 

7.24 ± 0.64 

8.48 ± 0.62 

3.27 ±0.23 

BridgeM 

7.70 ±0.11 

1.74 ±0.01 

7.19 ±0.39 

12.2 ±3.3 

3.35 ±0.79 

Bridges 

0.95 ± 0.04 

1.44 ± 0.04 

6.94 ± 0.90 

3.7 ± 1.1 

3.51 ±0.97 


“Parameters obtained by fitting a function 
F{E) = Ti(E/lGeV)“'"i exp{E/Ec) to Fermi-LAT data 
between 100 MeV and 300 GeV 
^’Parameters obtained by fitting a function 
F{E) = Fioo(F/ 100 GeV )“'"2 to MAGIG data between 

50 GeV and 400 GeV 



Figure 2: Light curve of the Crab pulsar at optical 
wavelength, 2.4 — 10 keV X-rays, 0.75 — 10 MeV, and 
100 — 300 MeV gamma rays (from top to bottom). The 
light curve at 50 — 400 GeV is overlaid on each plot for 
comparison. The optical light curve was obtained with 
the MA GIC telescope using the central pixel of the 
camera iLucarelli et al.ll20ok] . The keV and MeV light 
curves are from iKuioer et alJ [200 1|| . The 100 — 300 MeV 
light curve was produced using the Fermi-LAT data. All 
light curves are zero-suppressed by estimating the 
background using the events in the phase range from 0.52 
to 0.87. 

mechanism of pulsars. 
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Figure 4: Spectral energy distributions of the Crab 
nebula, PIm, P2m, Bridgeivi, and BridgeE measured with 
Fermi-LAT (below 50 GeV) and MAGIC (above 50 
GeV). The flux values averaged over the rotation period 
are plotted. 
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